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We establish the formalism in the nonsymmetric gravitational theory 



I/-) | (NGT) for stellar equilibrium and gravitational collapse. We study the col- 

lapse of a pressureless, spherically symmetric dust cloud. By assuming that 

^ the interior solution is smoothly matched at the surface of the star to the 

H , 

' quasi-static, spherically symmetric vacuum solution, we find that the star 

a' 

does not collapse to a black hole. It is anticipated that the finalcollapsed 

•i— I . 

| object will reach a state of equilibrium, and will emit thermal, gravitational 



and other forms of radiation, although the radiation may be emitted only in 
small amounts if the red shift from the surface of the compact object is large. 
No Hawking radiation is emitted and the information loss problem can be 
resolved at the classical level. 



From its seeming to me - or to everyone - to do so, it doesn't follow that it is so. What we 
can ask is whether it can make sense to doubt it. 

L. Wittgenstein, On Certainty 
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I. INTRODUCTION 



Recently, a new perturbatively consistent version of the nonsymmetric gravitational the- 
ory (NGT) has been formulated The linear approximation yields a theory free of 
ghost poles and tachyons and a Hamiltonian that is bounded from below |T|J2]. The flux of 
gravitational energy at asymptotic infinity is positive definite for a finite value of the range 
parameter a = /jT 1 . An analysis of spherically symmetric systems in the new NGT has been 
carried out by Clayton |J. 

It was conjectured that on the basis of a static spherically symmetric vacuum solution of 
the NGT field equations, in the long-range approximation, a 3> 2M, no black holes would 
form during the collapse of a star with a mass greater than the Chandrasekhar mass limit 
H|7|]. In order to establish the correctness of this conjecture, it is necessary to study the 
situation of the physical collapse of a star containing matter as the star's pressure falls to 
zero. 

In order to make the solving of the field equations manageable, we must make some 
simplifying approximations. We assume that the direct coupling to the skew part of the 
nonsymmetric source tensor is small and can be neglected during the collapse. We also make 
some simplifying approximations about the nature of the skew symmetric contributions to 
the field equations. With these approximations, we find that if the interior solution of the 
star and the exterior quasi-static solution are matched, the star collapses without forming 
an event horizon (black hole) and that the collapse can stop at some time, t = to, before it 
forms a state of singular energy density. 

The static vacuum (Wyman) solution contains no event horizons and the non- 
Riemannian geometry is singularity-free in physical spacetime. The NGT does not possess 
a Birkhoff theorem but it is reasonable to suppose that the final state of collapse 

will, under physically realistic conditions, tend asymptotically with time towards a pseudo- 
stationary final state. This state is one in which spacetime is invariant under an isometry 
group generated by a Killing vector field which is timelike for < r < oo. One anticipates 



2 



that nonstationary motions will in general be damped out by gravitational radiation, vis- 
cosity, etc. If such a stationary, massive final collapsed object forms, it will for practical 
purposes be indistinguishable from the observed massive "black holes" reported abundantly 
in the literature. The final collapsed object in NGT will radiate thermal radiation, gravi- 
tational radiation, etc., but not Hawking radiation, since the latter is a unique feature of a 
black hole event horizon. This could eliminate the problem of information loss associated 
with black holes in GR fH]. 

We begin, in Section II, with a presentation of the NGT field equations and the conser- 
vation laws. Then, in Section III, we provide a derivation of the motion of test particles, 
based on the fluid energy momentum tensor. In Section IV, we study the formalism for 
a spherically symmetric system and derive the time dependent field equations for such a 
system. In Sections V and VI, we develop the formalism for stellar equilibrium and gravi- 
tational collapse and in Section VII, we carry out an analysis of the field equations for the 
gravitational collapse of a spherically symmetric, pressureless dust cloud. This model of 
collapse is the equivalent, in NGT, of the Oppenheimer-Snyder collapse model in GR |3|. 
In Section VIII, we study the static spherically symmetric vacuum solution, and discuss 
the dual roles of the Riemannian and non-Riemannian geometries in NGT. In Section IX, 
we study the matching of the interior and the exterior metrics during collapse. Finally, in 
Section X, we summarize the results of this paper and discuss the claim made by Burko and 
Ori flH],|ll]] that black holes are anticipated in NGT collapse, using the linear approximation 



of an expansion of the NGT field equations about the Schwarzschild solution of GR. 

II. STRUCTURE OF THE NONSYMMETRIC GRAVITATIONAL THEORY 

The NGT is a geometric theory of gravity based on a nonsymmetric field structure with 
a nonsymmetric fundamental tensor g^ u , defined by 

9/iu = 9{fiu) + 9[fj,u]- (1) 
The affine connection coefficients, , are also nonsymmetric: 



— r^-j + r^. (2) 

We define the inverse tensor g^ w by the relation 

g" v g m = g Vil g«» = C (3) 

The NGT Ricci curvature tensor R^ U (W) is given by 

R,u{W) = - \{W^ U + w! Ptll ) - W^Wfr + W%W^ (4) 

where the are the unconstrained nonsymmetric connection coefficients, defined in terms 
of the affine connection coefficients through the relation: 

W£, = T% - |*JW„, (5) 

where Wp = W^ a] . It follows from (g) that 1^ = T^ x] = 0. The NGT Ricci scalar is given 
by R(W)=g^R^(W). 

The NGT field equations take the form: 



G^(W) + Ag^ + = 8vrT^, (6) 

2' 



rlH = _I g (^)w CT , (7) 



where 



and 



Moreover 



C^ v = -g^g [al3] g[0 a ] + g [a0i g^gp v + #[h> 



and P* = P^v — \g^vP with P^ = W^Wy. A and /i denote the cosmological constant and 
a "mass" parameter associated with the antisymmetric field g^u], respectively 
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We can write the field equations (0) in the form: 

R, U (W) = + 8tt(T^ - \g, v T - -^%„ + -^-P^), (8) 

Z OZ7T 4S7T 

and T = g^T^. Also, we have 

V = - ^g^C = g [af3] g» a gpv + ^g^g^^g 1 "® + g^u], (9) 

where C = g^C^. 

In empty space, the field equations (§) become: 

^V( r ) = + A ^ - + ( 10 ) 

where 

From the variational principle and the general covariance of the Lagrangian density, we 
can obtain the four Bianchi identities: 

\{G pv {T)z™ + G up {T)z™), a + l -G TU (T)g™ p = 0. (11) 

The Bianchi identities lead to the four conservation laws {§: 

^1*% + ^T"", + 2[/xi/, A]T^ = 0, (12) 

where 

[/^, ^] = g^M." + 9Xv,fi - g^,x)- 
This is known as the generalized law of energy-momentum conservation in NGT. 

III. EQUATIONS OF MOTION OF TEST PARTICLES 

Let us consider the equations of motion of test particles derived from the conservation 
law (0). We shall assume that the particle is confined to a tube E, whose linear cross section 

5 



dimensions are small compared to the length characterizing the gradient of the background 
metric. The fundamental tensor g^ v consists of a piece gj® corresponding to the continuous 
field at points along the world line of the test particle, and the part 5g^ u that describes the 
correction to the background gj® field due to the test particle |12j. We can then write 

9p.v = 9$ + dg^u- (13) 

To the first order of approximation, we keep in the field equations only the terms which are 
linear in 8g^ v . The energy-momentum tensor associated with the test particle is 8T^ U . It is 
assumed that the energy-momentum tensor T^ '^ u , associated with the background metric 
9uj}> vanishes inside as well as near the test particle. We adopt the notation: 

TV = y /=gST fU '. (14) 

The energy-momentum tensor for a fluid is 

= (p + p)vPu v - pg^ + K l " u] , (15) 

where = dx^/ds is the four- velocity of a fluid element, normalized so that 

g^uW = g { ^u v = 1. (16) 

In the test particle limit, p — > and we get 

The equations of motion of the test particle are then given by 

7-0") v _|_ J ^ V } T {al3) = F 
\a/3j 

where 

r = s^ a (g[Z] KmM m + Kmpv] [H , «] (0) ) , 

and we have used the inverse symmetric tensor s^ v defined by 
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(17) 
(18) 



to raise suffixes. Moreover, 



- j 



Let us define the proper mass of the test particle by the equation 

dt 



J d 3 xT 00 = u° J d?x^gp = ra t — , 



where po is the proper mass density , mt denotes the test particle mass and dr is the proper 
time along the world line. The equation of motion now takes the form: 

d ( ) (0) 

— (m^)+m t M M V = mi f, (20) 
dr [ ap J 

where we have defined 



We have 



From the condition: 



d du^ dm t 

— {m t u^) = m t — h —7—w 1 . 

dr dr dr 



dZ^u" = 1, 



we get 



1 



ito&^iOk = C^u\ = u/, ff = 0, 



where | denotes covariant differentiation with respect to the Christoffel symbol. We now 
have 



Multiplying (|20|) by we get 



dmt 



We now obtain the equation of motion of the test particle (we drop the (0) notation): 



dr 



+ 



(21) 



We see that, depending on the model chosen for K^ v \ there could be a contribution due to 
the non-conservation of the mass, caused by an exchange of energy with the skew field g^u]- 

According to Eq. fl2~T|) , test particles fall in an NGT gravitational field independently of 
their composition, so the weak equivalence principle is satisfied in the new version of NGT. 
However, the strong equivalence principle is not satisfied in the theory, because the non- 
gravitational laws of physics are not the same in different locally Minkowskian frames of 
reference; the skew part of the connection is a tensor which cannot be transformed 
away at a point, in contrast to the Christoffel connection j x \. 

If we assume that is small and can be neglected, we obtain the geodesic equation 

of motion for test particles Jl3 : 

(22) 



du^ | /i 
dr \ a(3 



u a vP 



0. 



IV. THE FIELD EQUATIONS FOR A SPHERICALLY SYMMETRIC SYSTEM 



For the case of a spherically symmetric field, the canonical form of g^ v in NGT is given 



by 



(-a w\ 

-p /sin^ 
-/sin^ -psm 2 9 



\-w 7/ 

where a, /3, 7 and w are functions of r and t. The tensor g^ v has the components: 



)1U 
















ur — c*7 

j_ /csce q 



\ 



/3 2 +/ 2 P 2 +P 

/CSCe /3CSC 2 6> 
>+/ 2 /3 2 +/ 2 







w— 07 











We shall assume that w = and only the gps] component of guuA is different from zero. It 
can be proved that only the static solution for gvm satisfies the physical, asymptotically flat 
boundary conditions p|.|I~lf . 
We have 

_,„2Wfl2_L f 2X11/2 



= sinfl[(a 7 - w 2 )((3 2 + f 2 )] 

The vector can be determined from: 

W„ = --JLs» p g^\ a . (23) 

For the spherically symmetric field with w = it follows from (p3|) that = 0. The field 
equations (0) and (|7|) for the spherically symmetric system now take the simpler form: 

G^(T) + Kg^ + -^ 2 C^ U = 8vrT^, (24) 
g M „ = 0. (25) 

Let us now reexpress the conservation laws in a form suitable for calculations in a spher- 
ically symmetric system. Using the Einstein notation, the compatibility equations are given 
by 

g^- = _ + gTM + ^ _ g^r? = 0, (26) 



where 



= a% + at;, 



and 



A^. a = A^ a + A p r^ p . 



We have 



Re[(G p _,_(r) g ^-) ;CT ] = i[(G^(r)g^-). ff + (G w (r) g ^-) ;CT ] = o. (27) 
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This can be written as 



Re[(G p _,_(r)g^-) ;(T ] = Re[(^(r)g-), CT - l -G TV {V){^Vl p + g"T£,)] = 0, (28) 
which leads to the result: 

Re[(G p _,_(r)g CT+ ^) ; .] = \{G pu {T)g? v ),« + \{G vp (V)fr)* + \^ru{T)g T \ p = 0. (29) 



This is the same as the Bianchi identities (|TlD - From (|24j) , we have for the spherically 
symmetric conservation laws: 

Bb[(T^u-S" +v ~U=0. (30) 
For a comoving coordinate system, we obtain 



u 



o _ ^ „,r 



v/7 
We define 



</ = u « = u * = 0. (31) 



= g^gauT ap , (32) 

which is a Hermitian symmetric tensor, T^j, = T vpi when g^,, is defined to be Hermitian 
symmetric: g pv = g up . From @ and ( |T5| ) we get 



T = p - 3p + <7 M 1^ = p - 3p + 2/*T, 

where we have defined: j^ 23 ] = KJ sin#. 

The NGT field equations in the presence of sources are given by 

^ 1 l r/ lA9 , „ 9n a' A' 7' / a' 7' 

i2 u (r) = — A - - \(A') 2 + 4B 2 } + + — - 

1U ' 2 8 [y J J 4a 2 7 V2a 2 7 

(l' V 9 /a\ d / 7 d 1 A 

~ V27J + <% V2^J + 27V27 ~ 2a~ + 2 J 

1 2 a/ 2 

= 47ra(p-p + 2fK), (33a) 
{2fB-(3A'\' t {2fB-(3A'^ 

3a 2 7 



i? 22 (r) = i? 33 (r)cosec^ = 1 + ( 2fB ~/ A ) ' + ( 2/ g~^ ) («'7 + y«) 
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^oo(r) 



B(fA> + 2JB) d(2fD-£A\ 2fD-0A f . 



4n[3(p-p + 2fK), 

-A--(A 2 + AD 2 ) + ^A + 



2"" 8 V ~~ y ' 47"" ' 2aV27 2a 



dt\2a 



A7 + 



2a 



7/ 2 



2aV2a 2 7 2 



4^ /5 2 + f 2 
4n 1 (p + 3p~2fK), 



R[io]{T) = 0, 

#(10) ( r ) = - 
= 0, 



1 1 a. 1 k \ liA 1 
-A' I A) -\ ' BD 

2 4 Va 2 / 4 7 2 



i? [23] (r) = sin# 



A '-?''( i+ ?f? 



= — An f sin 6{p — p). 



sin 9 



(33b) 



(33c) 
(33d) 

(33e) 



(33f) 



Here, prime denotes differentiation with respect to r, A — dA/dt, and we have used the 
notation: 



A = \n((3 2 + f 2 ), 

P 2 + f 2 ' 
Pf-fP 



B = 



D 



P 2 + f 2 



(34a) 
(34b) 

(34c) 
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V. STATIC EQUILIBRIUM EQUATIONS 



Let us consider the static hydrodynamic equations for stellar equilibrium. For a static 
spherically symmetric system g^ u ,p,p and are functions only of the radial coordinate 
r. For a fluid at rest, Eq.([H]) determines the velocity components. To simplify the field 
equations, let us assume that the direct coupling of gp3] to is small and can be neglected. 
We also adopt the long-range approximation /i«0 and we choose A = 0. The field equations 
are given by 



— A" UA'Y + 4B 2 H h — - 

2 8 LV ; 1 4a 27U 2 7 , 

iira^p — p), (35a) 



^22(r) = i?33(r)cosec = 1 + ^ — J + ^ — — J (a 7 + 7 a) 



+ + 2/35) 



4a 

= 4vr/3(p-p), (35b) 

= 47r 7 (p + 3p), (35c) 

%o](r) = 0, (35d) 

%o)(r) = 0, (35e) 

« P3l( r) - sH(^y + >- + < + $ - - h 

= -47r/ sin #(p-p). (35f) 



The conservation laws (|30|) give 



p' = -l(p + p)(ln 7 y. (36) 
We can further simplify the equations by adopting two approximation schemes: 

r 2 »/(r), (37) 

and 
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f(r) > r 2 , 



(38) 



where we have chosen j3{r) = r 2 . For the first approximation scheme, we obtain for the 
11/ar, 22/0,00/7 and [23]// field equations: 

~V^ + T^(- + -) + ^T + I = MP ~ P), (39a) 
2a7 4«7 V a 7 / a^r 

\-— 2 + 7^(---)+N = MP-P), (39b) 
r z ar z 2ra V a 7 / 



7" 7 



-^(- + ^) + ^ + L = 47r(p + 3p), (39c) 
2cry 4«7 V a 7 / r«7 

a'/' + jT + _ f + 4 =M ,-p). (39d) 



4a 2 / 4a7/ rcry 2a/ ra/ ra 2 

Here, we have defined 

a'P 8//' 8/ 2 3/' 2 a'ff ff" 

Ur) = ±- + -LL. - -J— J — + —±±- - LL- (40a) 

a 2 r 5 ar b ar b 2ar 4 2a 2 r 4 ar 4 

a 2 r b ar b a^r b 2a 2 r 4 ar 4 ar 5 2ar 4 2a7r 4 

Ur) = iUL - (40c) 
2a7r 4 a7r & 

Let us now consider the following combination: ll/2a + 22/(3 + OO/27: 

a' 1 1 

+ ^ - + G = Svrp, (41) 



ra 2 r 2 ar 2 



where 



Eq.(f4lD can be written: 



G(r) = -[I(r)+L(r)]+N(r) 



rY . „ , G s 2 



Integrating this equation gives 



1 - 2A4(r 

where 



1_8tt(p- -)r 2 . (42) 
a J on 
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M(r) = 4tt dr'r' 2 p(r'), 
Jo 



(44) 



and 



p(r) = p(r) 



G(r) 
8tt 



(45) 



From (|36|), we have 



7' 



2p' 



(46) 



7 P + P 

By using Eqs. (|39b|) , ([43|) and (|46|) , we obtain the extended Oppenheimer-Volkoff equation 

imp!: 



*M = -^(i + ^)(i 



where 



P(r) = 4vr 



p(r) + 



47rr 3 P(r)\/ 2.M(r)\- 1 
+ M(r) ' ~ 



(47) 



/N{r)-G{r) 



(48) 



V 47T 

Eq.fl4"7j) reduces in the non-relativistic limit to the Newtonian equation for stellar hy- 
drodstatic equilibrium. It is the same as the Oppenheimer-Volkoff equation, except that 
p(r) is replaced by p(r) in J\4(r) and p(r) is replaced by P(r) in the third factor on the 
right-hand side. 

By using (|47[) we get 



7(jJ r z 

The solution that has 7(00) = 1 is given by 



(49) 



-2 / — [M(r') +4vrr 3 P(r')] 



2A4(r') 



(50) 



These approximate results hold for a weak f(r) field and for moderately relativistic systems 
such as neutron stars. 

Consider now the second approximation scheme that holds for large f(r) with f(r) ^> r 2 . 
We get for the 11/a, 22/(3, OO/7 and [23]// field equations: 
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f 2 _ f_ df_ jW_ Y_ _ 5fV 

2a/ 2 ~~af + 2fa 2 + Aa^ 2 + A^a 2 ~ Ion ~ 2af 4 



/V a'f'r 4 
+ - h 



/ 3 

a r 



+ 



'^3 



/' if 



a' 



af 3 2a 2 f 3 a 2 f 2 af 3 af 2 
f" a'f 1 i f' 2 



r 2 afr 2a^f 2a 2 r af 2a 2 f ar 2 2a^r 2af 2 



Ya' 77' 



+ 



4a 2 7 ' 2«7/ ' 2a7 4a7 2 ' cry/ 2 2a r yf' 

4r 2 



47r(p-p), 
4vr(p -p), 
47r(p + 3p), 



a'/V 4 /'V /" a'f if 3fr 



>„3 



'„3 



+ 



+ 



7 r 



2a 2 / 3 af 2af Aa 2 f 4aryf af 3 cry/ 2 a/ 2 

Q,/ r 3 j?// r 4 y j?/ r 4 

_| 1_ I 1_ U. 

a 2 f 2 af 3 2ajf 

By forming the combination: 11 /a + OO/7, we get 



4?r(p -p). 



f/2 



P"^,4 



^ + 



+ 



2a/ 2 af 2a 2 f 2a f af 2a 2 f af af 



1 3 

ar 



+ 



77 /' 



+ 



8f'r 3 


8r 2 


a/ 3 


af 


~. 3 


/v 4 \ 


af 


2a/V 



) =87r(p-p). 



Using we find the hydrostatic equilibrium equation: 



P 



af(p + p) 



( /7 2 + 2/r 3_ /V 4) 



/2 



«7' 5f 2 r 4 f'r A 



af 2af 2a 2 f 2a f af 

'<-3 8^.2 ^3 

+ 87r(p-p) 



a'fr 8f'r" sr 



2a 2 f af 3 af a 2 f 



Near r = 0, we get 



P 



a/(p + p) 



f 



f f 



/2 



«7' 



+ 8vr(p-p) 



.a/ 2a/ 2 2a 2 / 
We adopt the equation of state for a gas of highly relativistic particles: 

1 



P=- 3 P- 



Using (|55f) in (|54|) gives 



(Aafp\(f f 
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a'f 16vr 



V /' Ac*/ 2a/ 2 2a 2 / 
In the limit of large densities, p — > 00, we get 



-P 



/64vr\/a/\ 2 
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This equation has the formal solution: 



p(r) = exp 



64tt f r , ,a(r')f(r')p(r') 



+ const. (58) 



f'(r') 

Thus, in the limit of large p and for ///' < 0, at the center of a dense collapsed object, 
only repulsive forces occur which increase with increasing density. For either a compact 
object of high density or a very massive object, we expect that an equilibrium state can be 
achieved. 



VI. SPHERICALLY SYMMETRIC COLLAPSE EQUATIONS FOR DUST 



If we adopt the approximation scheme leading to the geodesic equation, Eq . fl22|) , for 
falling test particles, then we can use a comoving coordinate system with the velocity com- 
ponents: 

U ° =1, u r = u = u 4> = 0> 

and the time dependent metric in normal Gaussian form: 

ds 2 = dt 2 - a{r, t)dr 2 - (3{r, t){d6 2 + sin 2 6d<j) 2 ). (59) 

In order to simplify the field equations, we must make several approximations. We shall 
assume that the star collapses as a pressureless dust with p = and that the direct coupling 
term in the source tensor, ifiw], is small and can be neglected during the collapse. As in 
the last section, we also set A = and take the approximation p w 0. The field equations 
now take the form: 

1 ,„ 1 ., ,,,, „ 9 . a' A' a a 2 aA 

2a 8a 1 1 Aa 2 2a Aa 2 4a 

1 1 / 2fB-(3A' V a' / 2fB-/3A' \ 

(3 4p V a J 8/3 V a 2 J 



i(^)-^ fA+2i3D) = i ^ (60b) 
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-±A--(A 2 + AD 2 )+ — 
2 8 V J Aa 



1 a 
2a 



2 .-. 47r ^' 



1 • , lA'a 1 A . ■ 1 

A + A A - = 0, 

2 4 a 8 2 

1 / fA' + 2PB \' l (fA' + 20B)a' IB^2JB^J3A)_ 
~Af\ a ) ~8 a 2 ] + 4 ^/ 

l (M + 2/?D)a + 1 }A + fA + + wb _ D _ pA = ^ 
8 af 4/ 4/ 

From the conservation law (^), we obtain within our approximation scheme: 



From this result, it follows that 



0_ 

df 



(pPVa) = 0. 



(60c) 
(60d) 

(60e) 



(61) 



There are two approximate regimes that we can adopt in order to further simplify the 
set of field equations. In the first one, it is assumed that 



/?M)»/(r,t), 



(62) 



while in the second one, we have 



/(r,t)»/?(r,t). 



(63) 



In the first approximation scheme using (|62"D, we obtain from (|60a|) -( |50^ ) the equations: 
-0" p' 2 a'p'^ ~ 2 



(3 2(3 2 2ap 



n < i - a(3 



+ -— - + W = Anp, 



2a Aa 2 2a/3 

1 _ 1 ffF_ _ a^_\ P_ aP_ = 

P a \2p AaPJ 2(3 Aap ^ 

a P a 2 p 2 „ A 



2a p Aa 2 2P 2 
pf P'P ap' 



~ + 



+ 



p 2P 2 2aP 



+ Z = 0, 



fP 



l + 



P 



12 



P 2 + l f'F 



2pf 2f 2ap 2 2P 2 2a 2 P aP P 2apf 

&/ di'P f" a'f 
Aaf 2aP 2af Aa 2 f 



Anp. 



(64a) 
(64b) 
(64c) 
(64d) 



(64e) 
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Here, we have defined 



TIT/ s a'(3'f 2 (3"f 2 a(3f 
W(r,t) = -tt4^ + 



2a 2 (3 3 a(3 3 2a 2 (3 3 

a&f _ 5^p aff clff _JT iff? 3/" , , 

2a(3 3 2a[3 4 2a(3 2 2a 2 (3 2 a(3 2 a(3 3 2a(3 2 ' 1 &j 

X ( t) - ^ P a,(3 ' f2 i P" f2 P p i i ^ 2/2 ff " 

lr ' } " 2a(3 3 2a 2 (3 3 a(3 3 (3 3 2a(3 2 [3 4 (3 2 a a(3 4 

offf f 2 If J'' 1 3///3 3/f/T 

2a 2 /? 2 2/3 2 /3 2 2a(3 2 2(3* 2a(3 3 ' 1 ) 

Y ^ t >- 2 /3 4 2(3 2 + (3 3 /3 2 ' (b5cj 



Z(r,t) 



(3'f 2 h(3(3'f 2 a(3'f 2 2/3 ff ff 3/7 
(3 3 2[3 A 2a[3 3 (3 3 (3 2 2j3 2 

2a(3 2 (3 3 



, lx J J *H J J 

+ WZTv + -53— (65d) 



In the second approximation scheme, we have for the 11/a, 22//?, 00, (01) and [23]// 
components: 



3 



2/a 2 2a/ 2 a/ 2a 4a 2 2a/ 2a 2 / 

5/' 2 /3 2 , f /3 2 a//? 2 d/3/3 3(3' 2 AfW a'(3'(3 

+ ~ T7t ~ 7T-7t. + —7^ + TT^Jt, = 4np, (66a) 



2a/ 4 a/ 3 2a/ 3 2a/ 2 a/ 2 2a/ 2 a/ 3 2a 2 / 2 

1 a/3 f 2 f af a'f (3 f" f(3 
(3 4a/3 2/ 2 + / + 2a/ + 2a 2 / 2(3 af + 2[3f 

(3" a'(3' (3'f f 2 
+ 2^-4^-2^/ + 2^ = 47rp ' (66b) 

2a + 4a 2 + 2/ 2 / 2/ 4 + / 3 / 2 2/ 2 + f 3 P ' 1 J 

/'/ f af 5f'fP 2 f{3 2 af(3 2 2f(3'(3 

2f 2 f 2a/ 2/ 4 f 3 2a/ 3 / 3 

3(3(3' 2(3(3 f (3'(3 a(3'(3 

+ — S H" + WZT,. = °. ( 66d ) 



2/2 /:•■ /2 2a/ 

3/gg /3 2 a//? 2 3/W /'V l' 2 «W 
2/3 + 2 / 2 2a/ 3 2a 2 / 3 a/ 4 2a/ 2 2a 2 / 2 

2/ 2a/ 2a 2 / 3 f 3 4a/ a/ 3 / 4 4a 2 / 2a/ 2 a/ 2 / 2 ^ 1 6J 
Here, we have used 
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B « , D « - . (67) 



VII. ANALYSIS OF DUST COLLAPSE 

We shall simplify our model for collapse even further and assume that the density p is 
independent of position. To analyze the field equations, we shall follow the procedures given 
by Weinberg, Landau and Lifshitz, and by Misner, Thorne and Wheeler P^TS|. Consider 
first the approximation scheme determined by the condition (|62[). It is assumed that a 
solution can be found by a separation of variables: 

a(r, t) = h{r)R 2 {t), (3{r, t) = r 2 S 2 (t). (68) 

From Eq. fl64d|) , we get 



| -§ = ^(r,t)r. (69) 
If, during the collapse, we assume that Z(r,t) ~ 0, then from (|69| ) we find that 

R{t) « S(t). 

Eqs. ( |64a| ) and ( |64b| ) now become 

+ R(t)R(t) + 2R 2 (t) + R 2 (t)W(r, t) = An R 2 (t) p{t) , (70) 



r/i 2 (r) 



+ ;/ + R(t)R{t) + 2R 2 {t) + R 2 {t)X{r, t) = ^R 2 {t)p{t). (71) 



r 2 r 2 h(r) 2rh 2 (r) 
The metric line-element takes the cosmological Friedmann-Robertson- Walker form: 



ds 2 = dt 2 - R 2 (t) 



h{r)dr 2 + r 2 (d6 2 + sin 



(72) 



Thus, we assume that during the collapse the metric is approximately isotropic and homo- 
geneous. If we assume that W(r,t) ~ X(r,t) ~ 0, then we find from Eqs. ([70]) and ( [71]) the 
GR solution: 
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where 



h'(r) 1 1 h'(r) 

2k = — = — - — — + 



rh 2 (r) r 2 r 2 h(r) 2rh 2 {r)' 
is a constant. 

We observe at this point that by means of a Killing vector analysis, it can be proved 
that a time dependent solution of NGT cannot describe an exact homogeneous and isotropic 
spacetime, unless is identically zero Let us expand the metric g( pv ) as 

SW) (r, t) = gfj u) (r, t) + 5g {liv) (r, t) , (74) 

where gW v s denotes the homogeneous and isotropic solution of g^ v ) and 5g^ v ) are small 
quantities which break the maximally symmetric solution with constant Riemannian curva- 
ture. 

Eqs.(|64cD and (170) can now be written: 



26(r) + R(t)R(t) + 2R 2 (t) + R 2 (t)W(r, t) = AnR 2 (t)p(t), (75) 

1 Ait 

- R(t)R(t) + -R 2 (t)Y(t) = Y R 2 (t)p(t), (76) 

where 

h'(r) 



26(r) 



r/i 2 (r) 

Eliminating .R by adding ( [75]) and fl76|), we get 



R 2 (t) = -b(r) + ^-p(r,t)R 2 (t), (77) 



where 



p{r,t) = p{t) - -?-H(r,t) (78) 

07T 



and 
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H(r,t) = ~[W(r,t) + ±Y(t)}. 
We shall normalize R(t) so that 

R(0) = 1, 

and we define 

p(r,t) = p(r,0)R- 3 (t). 

Let us assume that the fluid is at rest at t — 0, so that R(0) = 0. We have 

8n 

o(r) = yP( r ,0), 



where 



p(r,0)=p(0)-^-H{r,0). 

on 



We now obtain the equation: 



i?(t) = -6(r)+-|-p(r, 0)^(4), (79) 

which can be written as 

R 2 (t) = b(r)[R- 1 (t)-l}. (80) 

This has the same form as the corresponding equation in GR, except that b(r) has the 
additional contribution from the inhomogeneous, r dependent quantity, H(r,0), due to the 
skew fields. 

Eq.(|80|) has the parametric solution: 

ry + sinr? 

~hJV' (81) 

R= i(l + cosr]). (82) 

This solution reveals that R(t) = for 1] = ir and t — to where 
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7T 7T 



1/2 



2^b(~r) 2V87rp(r,0), 

Thus, as in GR, p(r, t) — > oo as — ► and the fluid sphere with initial density p(r, 0) > 
and p = will collapse from rest to a state with infinite proper density in the finite time 
to, provided that p(r,t) > 0. But p(r,t) need not be positive definite, for it contains second 
order skew curvature contributions. If, for R(t) near zero, we have p(r,t) < 0, then the 
collapse could be stopped even for the approximation f(r,t) <C f3(r,t). 

However, as we shall see in the following sections, the results for j3(r, t) ^> f(r, t) are only 
expected to hold at the initial stage of the collapse for r ^> 2M. As the collapse proceeds 
and the star becomes more dense, we should use the second aproximate regime for which 
the condition (|63"D holds. As before, we assume that a separable solution is possible: 



a{r, t) = q(r)R 2 (t), /3{r, t) = r 2 R 2 {t). (83) 

Then, equation ( |66cj ) becomes 

R f f 2 hf 2 r A R A fr A R 4 

— = -Ann - — + h 

R f 2/ 2 2/ 4 p 

8r 4 R 2 R 2 2r i R 3 R 8fr*R 3 R 

T 2 7 2 ~ + P ' ( } 

We can learn about the behavior of NGT by considering the behavior of fl84]) . Near 
R{t) = 0, Eq.(§§ becomes 

Let us assume that as the collapse approaches R(t) = 0: 

^-|-4vrp>0. (86) 

If for a contracting star we have R/R > 0, and since by definition R > 0, then R > and 
R(t) will not pass through zero during the collapse. Thus, the proper circumference and the 
proper three- volume of the collapsing star remain finite at the end of the collapse. Moreover, 
the proper density p(t) is finite for the final state of the star. In contrast, in GR we have 
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f = -fp(0)R- 3 , (87) 

from which it follows that the fluid sphere of initial density p(0) > and zero pressure must 
collapse from rest in a finite time to to a state of infinite proper energy density. 

VIII. EXTERIOR STATIC VACUUM SOLUTION 

The metric line element in NGT has the form 

ds 2 = m^ v) dx^dx u ', (88) 

where m^) can have three forms in the long-range limit ji ~ 0: 

l m^ u) = g (fU/) , (89a) 
2 m (^) = gO"), (89b) 

a m o«o = Hf^h-i _ ( 89c ) 

where s = Det^^)) and (7 = Det(^ /Ltl/ ). The metrics 1 mr l j lU ), 2 m^ v ' and 3 m^ u ^ were obtained 
from a study of the Cauchy evolution of field equations in Einstein's unified field theory by 
Maurer-Tison pp||2T|] . There are three light-cones in NGT, corresponding to the propagation 



of different zero mass modes. These three metrics can describe the causal Minkowskian light- 
cone structure of the spacetime in NGT, although it is expected that when \i 7^ 0, the light 
cone structure of the spacetime will be somewhat modified, since the massive gr^i field will 
not propagate causal information along the light cone. The three light cones degenerate 
to the standard single light cone of special relativity, when NGT is expanded about the 
Minkowski background metric, and also when it is expanded about a GR background metric. 
At the Schwarzschild radius, r = 2M, the three light cones are regrouped by interchanging 
their overlapping status . 



In Born-Infeld non-linear electrodynamics [22], there are two kinds of electric fields, one 



of which is point-like and singular at r = 0, while the other one is finite at r = 0; for the 
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latter the electric current density is spread out over space. Similarly, for the spherically 
symmetric NGT vacuum solution, the Riemannian geometries associated with the metrics 
1 m, 2 m and 3 m are singular at r = 0, while the non-Riemannian geometry determined by 
the fundamental tensor g^ v is finite at r = 0. The skew field g^] is like a "medium" which 
diffuses the spacetime metric and makes the non-Riemannian geometry non-singular. The 
ratio 

_ 

Rnr ' 

where Rr and Rnr denote the Riemannian and non-Riemannian scalar curvatures, behaves 
like a "dielectric constant" in spacetime. 

It is not meaningful to ask which geometry is the "correct" one, since the Riemannian 
and non-Riemannian geometries both play dual roles in the description of spacetime. We 
shall adopt for convenience the definition of the line element: 

ds 2 = g {liv) dx»dx h '. (90) 

We will assume that M<1//j. It can be shown that the only solution which yields an 
asymptotically Minkowskian spacetime has w(r) =0 |5j,p^|. 

In the case of the long-range approximation of NGT, corresponding to fi ~ in the field 
equations, the exterior static spherically symmetric solution has the form (Wyman |23f): 



7extW = e ". ( 91a ) 



M 2 e-\l + s 2 ) fduV 
aext[T) ~ (cosh(H - cos(^)) 2 \dr) ' [J ' 

Pext(r)=r 2 , (91c) 

2M 2 e _! ^[sinh(az/) sin(6z/) + s(l — cosh(azz) cos(6z/))] 

/ext( r ) = 7 T7 — \ 7TTX\2 ' ^ 91d ^ 

(cosh(az/) — cos(oz/)) z 



where 



VTT^+i j , /vTT 

a = \ and 



2 V 2 

M is identified with the mass and s is a dimensionless constant which is different for different 

bodies and is related to the strength of the coupling of matter to the skew field g^u]- 
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The function u(r) is determined by the relation: 

e i/ [cosh(az/) — cos(fo')] 2 2 = cosh(azz) cos(bu) — 1 + ssinh(az/) sin(tV). (92) 

Two coordinate systems can be used to analyze the Wyman solution, one of which is the 
standard spherically symmetric coordinates: x 1 = r,x 2 = 8, x 3 = <ft, x° = t. Another useful 
coordinate system uses the coordinates: x 1 = u, x 2 = 9, x 3 = (ft, x° = t, where 

M 2 e~ u {l + s 2 ) 

aext[u} ~ [cosh(H-cos(H] 2 ' [ } 

2A-f 2 [cosh(az/) cos(6z/) — 1 + s sinh(azz) sin(bu)} 
fjext{u) = e1cosh(H-cos(H] 2 ' 

and with 7 ex t(^) an d fexti^) gi ven as above. 

We must choose a particular branch of a solution of Eg. (|9"2"D . This choice is made by 
picking the branch that will yield the positive-mass Schwarzschild solution as a limit for 
r — > oo. This branch begins at v = and extends towards negative v. In such a coordinate 
system, the asymptotic weak-field region is at v = (r = oo), while the "origin" occurs at 
u defined by P(u ) = 0. The particular value of u depends on the value of s, and for s — 1, 
we find numerically that Pq ~ —5.1667. 

For r > M and \s\ < 1, the metric functions in conventional spherical coordinates are 
approximated by (for /x ^ 0) H0]: 

7ext( r ) ~ 1 - ( 95a ) 

Qfext(r) « (l - — J , (95b) 
sM 2 

/extW~-^-^(l + M r)- (95c) 

We see that for large r, the solutions for 7 ex ^ and a ex t have asymptotically the same form 
as the Schwarzschild vacuum solution in GR, and / ex t( r ) vanishes exponentially fast as 
r — > oo. 

Near r = we can develop expansions where r/M < 1 and < \s\ < 1. The leading 
terms are 
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4 7o (l + 0(s 2 )) / v~_\ 2 n ((r_\ 



/ext(r) = M 2 (^4 - K + S |s| + G(, 3 ) j + I * - ^ I ^ V 2 + 0(r 4 ), (98) 
where 

7o = exp(^ +0( S )J.... 

These solutions clearly illustrate t/ie non-analytic nature of the limit s — > m £/ie strong 
gravitational field regime < r < 2M P,|T|,[2~i| . Thus, the Wyman solution cannot be 
analytically continued to the Schwarzschild solution of GR for arbitrarily small values of the 
parameter s in the region < r < 2M. This means that the expansion about a fixed GR 
background, in the approximation of linear gwj, cannot be considered valid in the strong 
gravitational regime: < r < 2M. This fact will play an important role in our derivation of 
a solution to the collapse problem in NGT. 

To see that the non-analytic behavior in s of the Wyman solution for < r < 2M is not 
a coordinate dependent result, we can use the coordinate invariant norm: 



\s\ + s 2 n/8 + 0(s 



3\ 



where is the timelike Killing vector at spatial infinity: 

(7,0,0,0). 

Since 7 never vanishes throughout the spacetime, there are no event horizons and the solution 
is not analytic to the Schwarzschild solution, for < r < 2M as s — > 0, in any coordinate 
frame. The redshift is finite for s 7^ and the maximum redshift determined between r = 
and r = 00 is given by 

1 

Zmax — —= — 1- 

The singularity caused by the vanishing of a (r) at r = is a coordinate singularity, 
which can be removed by transforming to another coordinate frame of reference fgjTj • 
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The non-Riemannian geometry is non-singular in the range < r < oo, since all the 
nonsymmetric curvature tensor invariants are finite in this range of r. For example, from 
the non-Riemmanian curvature tensor: 

11 pup\ L ) — 1 1 f^p^v 1 au L fip ~r 1 ap L pv^ 

we find that the Kretschmann invariant: 

K = R x ^(T)R x , up (T), 

where 

R\pvp(T) = g\ a R a pvp(F) 

is finite. On the other hand, the curvature invariants formed from the Riemann-Christoffel 
curvature tensor B X p Vpi defined by 



- D pup 



J MU \, P i^p],u {av} {fip} + {ap} {{lis 



are singular like ~ M 4 /r 8 near r = P,|2"5|. The Riemann-Christoffel curvature tensors 
built out of the two metrics 2 m and 3 m will also have this singular type of behavior near 
r = 0. 

In the v coordinate representation, 7 ex t( z/ ) vanishes at the point v = — oo, corresponding 
to r* = (r 4 + J 2 ) 1 / 4 = 0, which is a point outside the physical causal spacetime. Moreover, 
there is a singularity in the non-Riemannian curvature invariants at the unphysical point 
v = — oo . 



Let us consider the exterior static spherically symmetric line element given by the metric 
m^u) = 9(pu)- 

ds 2 = 7ext(^ 2 " «ext(^ 2 - /W^ 2 + sin 2 9d4> 2 ). (99) 

We see that 7 ex t(^) does not vanish in the range: — oo < v < 0. Thus, all the points of 
spacetime are described by timelike Killing vector fields in the region — oo < v < 0; they 
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can be causally connected and all particles can be at rest. This corresponds to the fact 
that there are no trapped surfaces (event horizons) in this range of v. There are no event 
horizons in the r-coordinate frame in the region: < r < oo. 

Test particles which follow radial trajectories according to the path equation [|J: 

or test particles that follow radial trajectories according to the geodesic equation, Eq. (p2|) . 
are not stopped at r = but can continue through into the unphysical vacuum manifold for 



r < [13]. However, we shall see that in the collapse of physical bodies, it is unlikely that 



matter can collapse through the point r = and we conjecture that this is never possible. 

IX. MATCHING OF INTERIOR AND EXTERIOR SOLUTIONS 

We must consider now the matching of the metric outside the star with the one in the 
interior of the star. Birkhoff 's theorem does not hold for the spherically symmetric vacuum 
solution of NGT, which makes the matching of solutions a more difficult task to solve than 
in GR. It can be proved for /i ^ that no monopole radiation can escape to asymptotic 
infinity |J. 

The exterior metric outside the star can be expressed in terms of the coordinates, f, 8, <j), i 
in the form: 

ds 2 = C(f, t)dP - D(f, t)df 2 - 2E{f, t)dfdt - f 2 (d6 2 + sin 2 9d^ 2 ). (100) 
We remove E by defining the new time variable: 

dil = e(f,t)[C(f,t)dt- E(f,t)df\. (101) 



The line element now takes the form 161: 



ds 2 = 7ext( f ; tjdt 2 ~ «ext( f > t) df2 ~ f2 ( de2 + sin2 Qd<p 2 ), (102) 

where 
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Text ft*) = e' 2 (f,i)C-\f,i), (103) 
a ext ftt) =D(f J i) + C~ 1 (f,i)E 2 (f,i). (104) 

Let us first consider the approximation scheme in which Eq. ( |52j ) is valid. This approximation 
is expected to hold for the initial phase of the collapse when / is small. In order to match 
the solutions at the surface of the star, we must convert the interior solution with the metric 



(j72|) into "standard" coordinates. We shall assume that in the initial phase of the collapse, 



the interior metric is approximately determined by the solution Eg . (|73|) . We choose 

f = rR(t), 6 = 6, = 0, (105) 
and use an integrating factor which yields 

V k J JQ(r,t) (1 _ M) VI - RJ V ' 

where 

1 - fcr 2 \V2 



(1 fr r ^ \ 1/^ 



Here the constant vq is set egual to the radius of the star in comoving coordinates. We now 
get 



R/l - kr 2 \ l l 2 
QKl^kr 2 , 
kr 2 ^ _1 



7ft t) = 7,(^-2 , (107) 



«fti)=( v l-^J , (108) 
where Q and R are functions of f and i. At the radius of the star, we have 

f = f = r i2(t) (109) 

*= — ir-* L/ — 3Vh— 5 > ( n °) 



fc / m _ Vl - -R 

7(^ = 1-^. (HI) 

■^•M 1 -$)"'■ (n2) 
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The exterior and interior solutions match at r = roR(t) if we have 

l(ro,t) = 1ext( f o,t), (H3a) 
a(f ,t) = a ext (f ,t), (113b) 
f(r ,i) = fext(roJ). (113c) 

We shall now expand the exterior time dependent solution as 

7ext( f > *) = 7ext( f ) + 5 7ext(r, *), (H 4a ) 
aext( f >*) = «ext( f ) + Sa ext(r,t), (114b) 
W*S*) = /ext(^) + */«t(*S*)- (H4c) 

We assume that 57 ex ^,5a ex ^ and 5/ ex t are small quantities that can be neglected during 
the collapse (quasi-static approximation). The a ex t (f) , 7 e xt (^) anc ^ /ext(^) are determined 
by the unique static solution of the NGT vacuum field equations. 
For large values of the star's surface radius, r >> 2M, we have 

/ 2M\- 1 

a ext (r ) = (1 ~ —) , (115) 

7extM = l-— , (H6) 
?" o 

sM 2 

/ext( r o) = -^-exp(-r /a)(l + r /a). (117) 
If, as we have been assuming, r <C a, then we get from ( |117| ): 

/extfo) = ^. (118) 

Let us assume that s is related to the strength of the skew field coupling to matter by the 
expression: 

9 



(119) 



where g is a coupling constant, Nb is the baryon number of the star and (3 is some dimen- 
sionless constant. We expect that for a non-zero coupling of the skew field to matter there 
should exist a non- vanishing static exterior solution for /. 
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The interior and exterior solutions fit for large ro when r = roR(t) if: 



kr 2 = . (120) 

TO 

This yields for large r the expression for the total mass of the star: 

4<7T 

M = yP(0)r 3 . (121) 

For r R(t) — > 2M, we must switch to the approximation scheme based on the condition 
(|63|), since the quasi-static Schwarzschild solution fails to be a solution of the NGT field 
equations for r R < 2M |24j], provided we make the reasonable assumption that the limit 
to the static solution from small time dependent perturbations is smooth. We assume as 
before that the solution can be expressed as a separation of variables, as in Eq. fl83|) . Then 
the line element takes the normal Gaussian form: 

ds 2 = dt 2 - R 2 (t)[q(r)dr 2 + r 2 (d6 2 + sin 2 Odxj) 2 )}. (122) 

To express this line element in standard form, we must use the integration factor: 

V = J dRJ r (R ) ro), (123) 

where T is a function that is chosen to remove the cross-term drdt. The matching of the 
exterior and interior solutions is achieved by use of the conditions ( |113aj )-( |il8c| ). 



Because there are no trapped surfaces (event horizons) in the exterior Wyman solution 
in the range < f < oo, we have 

7ext (2M)^0, a ext (2M)<oo, / ext (2M) ^ 0. (124) 

From this we can deduce that the matching of the interior and the exterior solutions will 
not produce a black hole event horizon in the final stage of collapse. 

A light signal emitted in a radial direction at i will have df/di determined by ds 2 = 
where 

ds 2 = 7(f, t)dt 2 - a{f, t)df 2 - f 2 (d6 2 + sin 2 9d(j) 2 ). 
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The light signal will be detected at a distant point r at a time: 

F = t + 



rr' 

/ df 


[«ext( r )" 


1/2 


lr R(t) 


L7ext( f )- 





Because of the bounds in ( |124j ), we find that v < oo and the collapse to the Schwarzschild 
radius can occur in a finite time. 
In GR, we have 



F ~ -2krl In 



kr 2 , 

W) 



const, 



so that as tqR — > 2M = kr$, it takes an infinite amount of time for the star to collapse to 
the Schwarzschild radius. 
For large r , we have 

kr 2 



R(t) 



« 1 



and the redshift is given by 

z ~ ^(^|P) 1/2 ~ r Vk(l - R(t' - f')(R(t' - f) 

which is asymptotically the same as the result in GR. 

As ro — > 2M, the redshift is determined by the formula: 



1/2 



A — A dt ■ 
z = — ; = — - r R(t) 



dt 



a ext( r o R ) 



7ext( r o#) 



1/2 



(125) 



The red shift z seen by an observer is zero when the collapse is observed to begin, increases 
slowly during the collapse but satisfies the bound z < oo, because Eq.([L24]) holds for < 
r R(t) < oo. 

In GR, the redshift becomes for r Q R(t) —>■ 2M: 

krl \ ~ x ft' 



z~21 



R(t)J 



cxp 



2Arg7' 



(126) 



and the red shift becomes infinitely large as the collapsing star approaches the black hole 
event horizon. 
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Since a black hole event horizon does not form in the collapse of a star, in NGT, if we 
formulate the collapse problem with our matching conditions, the star is not cut off from 
the rest of the universe, and it can continue to emit all forms of radiation. The collapsed 
star is expected to form a dense static object that could be a strong source of gravitational 
radiation and other forms of radiation. 

X. CONCLUDING REMARKS 

We have been able to find an approximate treatment, in NGT, of the stellar equilibrium 
problem. We found that it is possible to realize a static equilibrium state for a massive 
compact object, due to the strong repulsive forces at the center of the object produced by 
the skew fields. We also formulated an approximate treatment of the gravitational collapse 
of a star considered as a spherically symmetric, pressureless dust sphere. By matching the 
exterior vacuum solution to the interior solution, we have shown that black holes are not 
expected to form during the collapse, since the exterior solution for small time dependent 
perturbations can approximate the static Wyman solution, which does not possess any 
trapped surfaces. The collapse can be stopped before the center of the star is reached. To 
completely solve the problem of the collapse of a star, in NGT, we must find a solution 
of the full NGT field equations, preferably for /x ^ and including a suitable source 
contribution. Such a solution is expected to be found only by using numerical methods to 
solve the field equations. 

It is important to stress that the collapse of a star, in NGT, must be stopped before 
R(t) = is reached, since the spacetime becomes unphysical for R(t) < 0. Because of 
the absence of black hole event horizons, such unphysical behavior would be "naked" and 
it would destroy the physical Cauchy data and make NGT a non-viable classical theory 
of gravitation. Of course, it is possible that we would have to discover a quantum theory 
of gravity to fully comprehend the collapse to small distances. But since such a theory is 
not presently available, it would be more satisfactory to be able to settle the issue of small 
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distance collapse within the classical regime of NGT. 

It has been conjectured by Burko and Ori [[K]] that black holes should be anticipated 
in the gravitational collapse of stars in NGT. Their analysis relied entirely on the use of 
an expansion of the NGT time dependent field equations about a GR background to first 
order in small g^u]- Specifically, this background was chosen to be the static spherically 
symmetric Schwarzschild solution. Of course, with this assumption a black hole is expected 
to form, because of the existence of a no-hair theorem for a skew symmetric potential coupled 



to a Schwarzschild metric [lljl , and the existence of trapped surfaces in the Schwarzschild 
solution. 

Since Burko and Ori demanded that / be small throughout the collapse of a star, 
they used a quasi-static exterior metric that can be closely approximated by the static 
Schwarzschild solution for small enough /. This metric cannot be a solution of the NGT 
vacuum field equations for < r < 2M [[24[] . The exact static spherically symmetric vacuum 



solution, in NGT, is the Wyman solution which does not have event horizons for < r < oo. 
It is not analytic to the Schwarzschild solution for arbitrarily small values of the parameter 
s in the range < r < 2M. From this we conclude that a quasi-static exterior metric is 
expected to be the Wyman metric plus a small time dependent part, which can be valid 
for suitably chosen initial value data with small /, before the onset of collapse. When the 
interior solution is matched to this exterior solution, black holes are not expected to form 
during collapse. This invalidates the claim made by Burko and Ori that black holes can be 
anticipated in the collapse of a star in NGT. 

It is incorrect to claim any definitive results for non-linear gravity theories, such as 
GR and NGT, on the basis of the linear approximation. One would not use the linear 
approximation in GR to solve the collapse problem, because if we expand the metric about 
Minkowski space, then for r = 2M, the perturbative expansion fails to be valid, since the 
metric perturbation is of order 1. Similarly, the Burko-Ori quasi-static expansion of g^ v 
about the Schwarzschild background fails to hold at r = 2M, because becomes larger 
than unity and the linear approximation breaks down. Given that a generic source exists 

34 



that generates a non-zero g^ , a static solution for / is produced. Then, we know from the 
exact Wyman solution that / becomes large near r = 2M, and the linear equation for / fails 
to describe correctly the physical collapse of a star. Assuming that / is just a radiating wave 
without a static part during the collapse is not a physically realistic treatment of the problem 
in NGT. Burko and Ori incorrectly assumed that by adding matter and a generic coupling 
of matter to the skew fields for gravitational collapse would not change their conclusions. 

The NGT violates the strong equivalence principle and, therefore, it is expected that 
a freely falling observer (the weak equivalence principle is not violated in the new massive 
version of NGT) could perform experiments to detect the formation of an event horizon. 
This is not possible in GR. This effect would only show itself in higher orders of g^ , and 
would have important consequences for the physics of collapse in NGT. 

Actual stars would collapse more slowly than in the model which we have studied because 
of the effect of the pressure of radiation, of matter and of rotation. 

Since the final collapsed object is expected to be a massive compact star without an 
event horizon, radiation of all forms can be emitted by the surface of the star. Of course, 
if the red shift emitted from the surface of the star is too large (but never infinite), then 
in practice only small amounts of thermal and gravitational radiation will escape. There 
would be no Hawking radiation emitted, for such radiation is associated specifically with 
a rigorous black hole event horizon. Therefore, the problem of information loss associated 
with the quantum mechanical aspects of a black hole would be eliminated. 
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APPENDIX A: THE TIME DEPENDENT T-CONNECTIONS 



The NGT compatibility equation is given by 

g\u, v - gpu^xr, - g\pT^ = ^g M (g P ugx v - g w g\ P - g\vg\ m \W^ 

where is determined from fl23|). For the spherically symmetric system, when 
g[oi]{r,t) = 0, it follows that = and the compatibility equation reads: 

The non-vanishing components of the T-connections are: 

pi - a ' 
r i _ o 

T\ 2 = T\, C ose C 2 6 = ±-(fB - )^f3A 
ioo_ 2a' 

1 (12) - 1 (13) - 4 A > 

r 2 _ F 3 _ 1 4 
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r [is] = - r [i2] sin 2 9 = -B sin 0, 
r [so] = -rf 20] sin 2 9 = -\p sin 0, 
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where A, B and D are given by Eqs.( |34a|) , (|34b|) and ( |34c|) . 
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